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ABSTRACT 


For  high  frequency  monochromatic  waves  the  geometrical  theory 
of  diffraction  can  be  used  for  predicting  radar  scattering.  For 
typical  targets  the  strongest  returns,  due  to  speculars,  occur  only 
at  special  aspect  angles.  At  most  aspect  angles  the  dominant  returns 
are  usually  due  to  diffraction  by  the  "edges"  of  the  target. 

The  proposed  use  of  short  pulse  radars  requires  the  consideration 
of  scattering  by  non-monochromatic  signals.  The  progressing  wave 
formalism  is  a  generalization  of  the  geometrical  theory  of  diffraction 
suitable  for  treating  pulses  for  which  the  high  frequencies  predominate , 
This  paper  extends  the  earlier  work  on  specular  scattering  of 
progressing  waves  to  include  scattering  by  edges.  For  a  moving 
target,  partial  results  are  obtained. 
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SECTION  I 


INTRODUCTION 

The  "geometrical  theory  of  diffraction"  ^  is  useful  for  the 

solution  of  problems  of  scattering  of  time-harmonic  waves  at  high 

frequencies.  A  recent  generalization  of  that  theory,  the  "progressing 

[2] 

wave  formalism"  is  suitable  for  scattering  of  pulses  of  arbitrary 

form  by  fixed  or  moving  objects.  In  both  theories  the  largest  returns 
(specular  returns)  are  due  to  reflected  waves,  but  for  typical  objects 
specular  returns  occur  only  at  special  aspect  angles,  and  «tt  other 
angles  the  main  returns  are  usually  due  to  diffraction  by  edges. 

The  purpose  of  this  paper  is  to  derive  the  formulae  for  the 
diffracted  progressing  wave  produced  by  an  arbitrary  edge.  The 
reflected  wave  is  discussed  in  [2].  The  theory  is  based  on  the 
scalar  wave  equation,  but  the  generalizat  ion  to  the  electromagnetic 
case  should  be  straight-forward  (see  [3j  and  [A]) .  In  Section  II 
we  discuss  the  progressing  wave  formalism  for  the  wave  equation. 

This  discussion  is  a  simplified  and  improved  versicn  of  the  basic 
theory  presented  in  [2] #  The  use  of  analytical  formulas  (involving 
jacobians)  rather  than  geometrical  formulas  (involving  area  ratios) 
for  the  amplitude  coefficients  is  a  great  convenience.  An  improvement 
in  the  theory  is  also  achieved  by  the  use  of  finite  part  integrals, 
which  were  recently  exploited.  In  Section  III  we  derive  the 
formulas  for  the  diffracted  wave.  The  method  used  relies  heavily 
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on  the  geometrical  theory  of  diffraction,  but  there  are  several  new 
features  that  do  not  appear  in  the  time-harmonic  case.  A  special 
case  of  the  formula  derived  was  obtained  earlier  in  [5],  but  there 
only  "time-independent"  progressing  waves  (see  Section  II)  were 
treated  and  the  discussion  was  restricted  to  the  two-dimensional  case. 
In  Section  IV  we  obtain  partial  results  for  diffraction  by  a  moving 
edge.  (A  complete  theory  would  require  the  solution  of  a  non-trivial 
"canonical  problem") .  These  results  show  some  interesting  new  feature 
involving  the  direction  of  the  diffracted  ray  and  a  "Doppler  shift" 
in  the  instantaneous  frequency.  The  appendix  contains  a  discussion 
of  the  fractional  integration  operator,  which  is  essential  to  the 
theory  of  edge  diffraction. 
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SECTION  II 


THE  PROGRESSING  WAVE  FORMALISM 


Let 


,  v  -iOJt  ,  . 

v(t  ,  x)  =  e  u(x) 


(1) 


be  a  solution  of  the  wave  equation 


Lv  =  v 


xx  2  ,  .  tt 

V  V  c  (x) 


=  0 


(2) 


Here  uu  is  a  constant  and  the  repeated  index  V  is  summed  from  1  to  3, 
By  inserting  (1)  in  (2)  we  find  that  u  satisfies  the  reduced  wave 
equation  (for  simplicity  we  take  c  to  be  constant) 


u  +  k2u  =  0 ,  k  =  - 


(3) 


Contemporary  generalizations  of  geometrical  optics  are  often  based 
on  asymptotic  solutions  of  (3)  of  the  form 


,  iks  (x)  \  , . ,  >.  -m  (  v  , 

u(x)^e  —  /  (ik)  z  (x)  ,  k  “• 

—  La  m  ~ 

m=0 


(A) 
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If  we  insert  (4)  into  (l),  multiply  by  an  arbitrary  function  af'uu), 
and  integrate  with  respect  to  0),  we  obtain  a  formal  solution  of  (2)  of 
the  form 


v(t  ,  x)  =  )  e  [t  -  (x)  .  (5) 

~  <_i  m  c  — 

m=0 

Here 

z^Qc)  =  (-c)mzm(x)  (6) 

and 

e  (t)  =  I  a(uu)(-iuu)  me  ^^dU)  .  (7) 

m  v 

We  note  that  the  functions  e  (t)  are  successive  anti-derivatives,  i.e., 

m 

e  (t)  =  e  (t)  .  (8) 

m  m- 1 

A  solution  of  (2)  of  the  form  (5)  is  called  a  time- independent 
progressing  wave,  provided  (8)  is  satisfied.  Since  (4)  is  valid  for 
large  k  *  U)/c,  it  is  clear  that  (5)  is  valid  provided  the  spectrum 
a(uu)  contains  only  high  frequencies,  or  at  least  a  (ou)  is  small  except 
for  large  UU.  This  provides  an  interpretation  of  the  formal  expansions 
we  shall  construct.  Other  interpretations  are  discussed  in  [2]  , 

We  now  generalize  (5).  A  (general)  progressing  wave  solution 
of  (2)  is  given  by  an  expansion  of  the  form 
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00 


v(t,  x)  “  )  e  [0(t,  x)]  z(m)(t,  x)  (9) 

m=0 


provided  (8)  is  satisfied.  The  function  eQ(t)  is  called  the  wave- form 

and  is  essentially  arbitary.  We  assume  that  all  of  the  e  (t)  vanish 

m 

at  t  =  -co.  Then  the  successive  e^(t)  are  determined  uniquely  by  (8)  , 
In  order  to  determine  the  phase  function  0(t,  x)  and  the 
amplitude  functions  z^m^(t,  x) >  we  insert  (9)  into  (2)  and  use  (8), 
This  yields  the  characteristic  equation 


0  0 
xv  \ 


0 


(10) 


and  the  transport  equations 


2  [0 


z 

X  X 

v  v 


(m) 


1  m  (®)l  W,*  T  (m-1) 

—  ^tzt  J  +  2  10  =  -  Lz  , 


m 


0,  1,  2 


0  . 


(ID 


We  introduce  the  functions 


-  -  0fc(t,  x)  .  \  =  0X  (fc,  i)  (12) 

v 


5 


and  set 


k  *  (ki ,  k2  »  ^3)  ,  k  *  Vkv\> 


03) 


Then  (10)  takes  the  form 


2  2,2 

uu  =  c  k 


(14) 


or 


W  *  hfk,  x)  *  :fc  c(x)k 


05) 


The  first  order  partial  differential  equation  (15)  can  be 
solved  by  che  method  of  characteristics.  We  introduce  the  character¬ 
istic  equations  (Hamilton's  equation) 


x  =  g 
V  6V 


dh 


(V  -  1,  2  3) 


(ib; 


where  £  =  (g^,  ,  g^)  is  the  group  velocity  vector  given  by 

gy  =  |~-  =  ±  ck^/k  *  c^/U)  (v  *  1,  2,  3)  .  (17) 

A  solution  x(t) ,  k(t)  of  the  system  of  ordinary  differential  equations 
(lb)  determines  a  curve  x  *  x(t)  in  x-space  called  a  ray .  This  ra  ' 
is  the  projection  onto  x-space  of  the  space-time  curve  [t,  x(t)3. 

The  latter  is  called  a  characteristic  curve  of  (10)  ,  or  bi character is tic 
of  (^2)  .  We  note  that 
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(18) 


ub  «  dh/dk^  +  dh/dx^  x^  =  0 


and 


i.e.  U)  and  0  are  constant  on  bicharacteristics.  Thus  the  bicharacter¬ 
istics  lie  on  level  surfaces  of  0(t,  x)  .  These  level  surfaces  are 
called  characteristic  hyper surfaces  of  (2) .  If  the  value  of  0  is 
prescribed  on  some  hypersurface  or  lower  dimensional  manifold  M,  the 
values  of  U),  and  k  can  then  be  obtained  by  differentiation  along  M. 

Then  the  bicharacteristics  emanating  from  M  can  be  obtained  by  solving 
the  initial  value  problem  for  (16) 5  and  from  (19)  the  value  of  0  on 
each  bicharacteristic  is  equal  to  its  value  on  M.  This  procedure  for 
determining  the  value  of  0(t,  x)  will  be  illustrated  in  Sections  III 
and  IV.  It  is  also  discussed  in  detail  in  [2]. 

In  order  to  determine  the  function  z*m^  we  write  (11)  in  the 


form 


V 


(20) 


Thus  we  see  that  the  transport  equations  are  in  fact  first  order  linear 
ordinary  differential  equations  along  the  bicharacteristics.  A  brief 
digression  is  now  required  in  order  to  obtain  the  solutions 
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in  a  convenient  form.  Let  us  assume  that  we  have  a  3-parameter 
family  of  bicharacteristics  [t,  x(t ,  a)]  where  a  ■  (a^,  ,  O is 

the  parameter  vector  (e.g.  a  might  be  the  parameter  vector  in  a  parametric 
representation  of  the  hypersurface  M)  .  We  introduce  the  jacobian 


j  -  j(t,  c) 


-  [  £  U.  2>  ] 


d(xis  x2,  x3) 

3(a1,  a2,  o3) 


(21) 


and  observe  that 


z 


V=  1 


da 


cof 


V 


(22) 


Here  6^  is  the  Kronecker  symbol  and  "cofM  denotes  the  cofactor.  If 
i  =  k  (22)  follows  from  the  rule  for  the  expansion  of  a  determinant 
by  cofactors.  If  i  ^  k,  it  follows  from  the  fact  that  a  determinant 
with  two  identical  rows  vanishes.  By  differentiating  the  determinant 
we  find  that 


(23) 

Here  we  have  used  (22)  and  (16).  We  now  set  f  ■  (f^,  f ^ ,  f^)  vher* 

f.  -  k./<«  -  g./c2  ,  i  -  1,  2,  3  .  (24) 

it  i 
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Since  UJ  *  -  0 

*1  "1 


(k^)  ,  Vo)  =  -  k,_  ,  hence 


19  12  2  -2 

kt.k  =  J  (k  ) t  -  j  (ur/c  )t  -  c  uu>t 


(25) 


Therefore 


;  [7-i +  7  “t]  -  5  cVv-  7  4«]  -  i  w  • 


(26) 


It  follows  from  (23)  that 


3J  log  hi  -  if*  (<=2f>  -=2  v.f+£.^ 


5T  w  +  T  ^ 

c 


(27) 


We  now  see  from  (16)  that 


k1**  ill 


l_  d_ 
2  dt 


2 

c 


c 


d_ 

dt 


log  y 


2  dy_ 
y  dt  5 


(28) 


where 

y  -  |j|*/c  .  (29) 
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By  inserting  (28)  into  (20)  we  obtain 


This  equation  can  now  be  integrated  to  yield  z^m^(t)  in  the  form 


y(t)z^m\t)  =  y(t1)z^(t1) 


1  f  2  _  (ra-1),  . 

2i  J  c  yLz  dt 


(31) 


Here  we  have  not  indicated  explicitly  the  dependence  on  the  parameters 

C,  Thus,  e.g,  z(m)(t)  =  z(ra)  [t,  x(t ,  0_ )  ]  is  the  value  at  time  t  of 

z(m)  on  the  bicharacteristic  Ct,  x(t,  a)]. 

If  j(t)  vanishes  at  t  =  tQ,  the  space-time  point  [tQ,  *  (tQ,  0)3 

is  called  a  caustic  point.  Then  y(t  )  vanishes  and  z^(t  )  becomes 
- * -  o  o 

infinite.  If  t  =  t^  is  a  caustic  point  (31)  is  not  valid  because 
z(tn)(tl)  is  infinite  and,  in  general,  the  integral  diverges.  In 
order  to  obtain  a  valid  formula  for  this  case,  we  first  define  the 
"finite  part"  of  a  function.  For  t>t  let  f(t)  have  an  asymptotic 
expansion  in  powers  (perhaps  fractional)  of  (t-tQ)  as  t  ■*  tQ.  Let 
f  (t)  denote  the  singular  terms  [negative  powers  of  (t-t  )  ]  in  this 
expansion.  We  define  the  finite  part  of  f(t)  as  t  •*  t  by 

fin  f(t)  *  lim  [f(t'>  -  fw(t)]  (32) 
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a 

J 


Now  if  J  g(x)dx  is  divergent  or  convergent  at  x  ■  t  we  define  the 


finite  part  of  the  integral  as 


g(x) dx  «  fin  J  g(x)dx  . 


t  t 


o  t 


(33) 


We  now  set 


(m) 


(m) 


Cv  =  fin  y(t)zv  ' (t)  =  fin 


t  -»  t 


t  -  t 


luolV^u) 

c(t) 


(34) 


and  take  the  finite  part  as  t^  “*  t  of  (31)  .  This  yields 


y(t)z(m)(t)  =  CU)  -  ~  ^c2yLz(m_1)dt-  , 


(35) 


or 


■  (m)(t)  =  ^(m)c^)  -  ^  <pc(t)c(tl)  I'i^-|\z(m_1)(t,)dt'  , 

h(or 


m  =  0,  1,  2,  z(_1)  =  0  . 


(36) 


If  t  *  t  is  not  a  caustic  point  j(t)  and  z^m^(t)  are  continuous 
there.  Then  (34)  yields  =  y  (tQ)  z  ^  (t  )  and  the  finite  part 

integral  becomes  an  ordinary  integral.  In  this  case  (35)  or  (36) 


reduces  to  (31)  with  t. ^  replaced  by  t  . 
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In  (36)  c(t)  -  c[x  (t ,  a)].  For  the  special  case  of  a  homogeneous 
medium  c  is  constant  and  (36)  becomes 


>>(t) 


.(m) 


g(t) 


c2  L|J(t')i  (m-1),  ,.v.  , 

T  '  55  <  I  rn  ' 1  L  (t  )dt  ' 


m  “  0  ,  1  J  2  y  «  «  »  y  z 


(-D  . 


Here 


5(">  -  fin  b(t)|V"\t>  . 

t  t 


(37) 


(38) 


When  c  is  constant  we  see  from  (16)  that  k  is  constant  on  a  bicharacter¬ 
istic  and  from  (17)  that  £  has  the  constant  value 

g  =  c  u 


where  u  is  the  unit  vector 


(40) 


Then  the  first  equation  (16)  can  be  integrated  to  yield  the  explicit 
ray  formula 


x  *  x(t  i  2)  “  x  (2.)  +  cCt-t  (a)  ]u(a)  . 


(41) 


The  bicharacteristic  [t,  x(t ,  a)]  is  now  a  straight  line  in  space 
time.  It  passes  through  the  point  [t  (a),  x  (a)]  . 
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SECTION  III 


DIFFRACTION  BY  EDGES 

Let 

v  (t,  x)  =  Y  e  (0)z^  (42) 

o  —  Lj  m 

m=0 

be  a  progressing  wave  solution  of  (2)  with  constant  c.  Let  S  be  a 
boundary  surface  in  x-space ,  such  as  the  surface  of  a  fixed  scattering 
obstacle.  In  space-time  the  boundary  is  a  hyper-cylinder  S  with 
generators  parallel  to  the  t-axis.  The  projection  of  S  on  x-space 

/N/ 

is  S.  If  the  bicharacteristics  of  (42)  intersects,  a  reflected 
progressing  wave  is  produced.  Such  problems  are  discussed  in  detail 
in  [2].  In  an  important  special  case  the  incident  wave  is  a  plane 
wave 

vo(t,  x)  -  e0(t  -  x1/c>  <43) 

where  e^  is  an  arbitrary  function.  Then  (43)  is  not  only  a  special 
case  of  (42),  but  is  in  fact  an  exact  solution  of  (2), 

An  edge  E  is  a  curve  along  which  S  fails  to  be  smooth.  Examples 
are  the  edge  of  a  screen  or  the  edge  of  a  finite  cylinder  or  other 
object  which  is  locally  wedge-shaped.  If  S  contains  an  edge,  we 
assume  that  in  addition  to  the  incident  and  reflected  waves  the 
solution  contains  a  third  term,  the  diffracted  progressing  wave 
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A  r1  A  A  A 

v(t,  x)  -  )  «  (0)zW  .  (44) 

m*0 

A 

We  assume  that  the  phase  function  0  of  the  diffracted  wave  is  equal 
to  the  phase  function  0  of  the  incident  wave  on  E,  Let  r|  be  an  arc- 
length  parameter  and  x  *  x  On)  the  parametric  equation  of  E.  Then 


0[t,  XjOl)] 


(45) 


By  differentiating  (45)  with  respect  to  t  and  T)  we  obtain 


and 


A 

U)  *  OJ  on  E  , 


(46) 


A 

k.t^  •  k.  t^  on  E 


(47) 


A  A  A  A 

Here  UU  «  -  0^ ,  k  ■  V0 ,  and  t^  *  dx^/dr]  is  the  unit  tangent  vector  to 

q  A  C  ^ 

E.  In  terms  of  the  unit  vectors  u  •  —  k  and  u  »  r-  k,  which  point  in 

—  (JD  *—  *—  m  ' 

CL) 

the  directions  of  the  incident  and  diffracted  rays,  (46)  and  (47)  yield 


u«  t^ 


cos  0 


(48) 


*  where  (3  •  P(t,  T))  is  the  angle  between  the  incident  ray  and  the 


edge,  i.e. 


cos  (8  •  — 1 


(49) 
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This  is  the  familiar  law  of  edge  diffraction  (see  [ l] ) •  If  we 
introduce  the  unit  normal  and  binormal  vectors  n  and  b  of  E  we  see 
from  (48)  that 

A  A 

u  =  u(t,  T),  a)  =  cos  @t^  +  sin  (B  cos  an  +  sin  (3  sin  a  b  . 

(50) 

This  equation  defines  the  angle  a.  In  order  to  satisfy  (45)  we  find 
that  at  each  time  t  =  T  a  one  parameter  family  of  diffracted  bi¬ 
characteristics  emanates  from  each  point  XqCt))  of  E#  The  corresponding 
rays  are  given  by 

r  n  A 

X  =  x(t ,  T ,  tj,  a)  =  +  c[t  -  T-J  u(T,  a)  >  0^a<2rr  %  (51) 

They  generate  a  cone  with  vertex  at  and  semiangle  3.  Equation  (51) 
is  a  special  case  of  (41)  with  parameters 

(a1,  a2,  a3)  =  (t,  ri,  a)  .  (52) 


From  (18),  (19),  and  (40)  we  find  that  on  the  bicharacteristic 
(51), 

0)  *  uu[t  ,  x  On)  ]  ( 53) 

— o 

0  =  0[t  ,  x  f Tj )  3  (54) 
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and 


A  yj  A 

k  -  S  u(t,  ri ,  a)  .  (55) 

In  order  to  use  (37)  to  determine  the  amplitude  coefficients 
we  must  compute  the  jacobian 

j(t)  -  j(t,  T,  n,  a)  -  d(xl>  *2>  X3}  (56) 

3  (t  ,  ti  ,  a) 

from  (51),  In  so  doing  we  make  use  of  the  Frenet  equations, 

t  1  =  Hn  n*  *  -  ^t  +  T  b  ,  b1  «  ■  T  n  ;  1  «  ^  .  (57) 

i  -  1  -  — 1  o-  *  —  o-  dr) 

Here  X  and  T  are  the  curvature  and  torsion  of  the  curve  E#  A 
o 

straight-forward  computation  yields 

j(t)  -  cosin2?  (1  +  |)  -  c2sin2e  (t  -  T)  [l  +  '  T^],  (58) 

where 

a-c(t-T)  (59) 


is  che  distance  along  a  diffracted  ray  from  the  edge. 


sin2p 


3'sinP  +  —  sinf3cos£  +  KCos6 
c 


,  3'  *  3^(7,  H)  ,  3  -  3t  (T,  TT)  , 

(60) 
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and 


cos&  =  sin|3cos&  =  u.n 


(61) 


From  (37)  with  m  =  0,  and  from  (58)  we  now  see  that 


Z(0)(t)  =  y(0)(t,  t|,  a)  |o(l  +  2)  |  ,  a  =  c(t  -  t) 


(62) 


Here  we  have  set 


,(0) 


-(0) 


c  2sin0 


(63) 


Thus  in  order  to  obtain  the  leading  term, 


A  A  A  a  /q\ 

v(t ,  x)  ~  e  (0) 


(64) 


of  the  diffracted  wave  it  remains  to  determine  the  functions  e  (t) 

o 

(0)  A (0) 

and  Y  (t,  T],  Ol)  ,  We  assume  that  the  amplitude  z  of  the  diffracted 

wave  is  proportional  to  the  amplitude  z^^  of  the  incident  wave  at 

the  (space- time)  point  of  diffraction.  Thus  we  set 


y(0)(t,  n,  a)  =  d  x  On)]  . 


(0), 


— O 


(65) 


Here  the  diffraction  coefficient  D  is  to  be  determined. 
-  o 
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The  value  of  depends  on  the  boundary  condition  at  S  and  on 
the  local  geometry  at  the  edge.  We  consider  the  two  boundary 


conditions , 


(A)  v  =  0  , 


on  S 


(66) 


illustrated  in  Figure  1.  The  plane  of  the  figure  is  orthogonal  to 
the  edge  and  the  unit  tangent  vector  t^  points  into  the  plane.  The 
projections  of  the  incident  and  diffracted  rays  onto  the  plane  of  the 
figure  are  shown,  y  is  the  local  wedge  angle  which  is  bounded  by 
the  tangents  to  S  at  E.  The  angles  9  and  a  are  simply  related  as 
shown  in  the  figure,  i.e. 


9  =  Y  -  a 


(67) 


where  Y  =  Y(r|)  is  the  angle  between  the  unit  normal  vector  n  and  the 


vector  t^  which  is  orthogonal  to  one  tangent  plane  of  S  at  x^ri). 

A 

In  order  to  compute  e^  and  Dq  we  begin  with  an  incident  wave 


of  the  form  given  by  (1)  and  (4).  Then 


00 


i[ks  (x)  -UUt] 
v  ~  e  — 


k  =  uu/c  , 


(68) 


o 


m=0 
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Figure  I .  Rays  of  Edge  Diffraction 
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and  according  to  the  geometrical  theory  of  diffracction 
the  diffracted  wave  is  given  by 


[i] 


A 

V 


00 


-h  itks  (x)  -U)t ] 
e  — 


(x)  . 


(69) 


The  leading  term  is  given  by 


where 


A 

z 

o 


ji^/^Dz  [x  (r|)]  |a(l  +  —) 
o  o  p 


(70) 


sin(TT/  q) 

(2TT)  ^qsin(3 


e-9o  I -V  |  n 

cos  -  I  +  Icos  — 

q  1  \  q 


cos 


0+0  +TT 
o 

q 


]• 


(71) 


v 

Here  q  =  2  -  —  and  p  is  given  by  (60)  with  (3  *  0.  The  upper  or 
lower  sign  in  (71)  holds  for  the  boundary  condition  (66)  (A)  or  (B)  # 
We  now  set 


0)  =  bV  , 


A  A 

0  a  b(t  -  s/c)  ,  0  *  b(t  -  s/c) 

(72) 


where  b  is  a  constant.  Then  k  *  bv/c  and  we  may  re-write  (68)  and 
(69)  in  the  form 


v 

o 


X 


-iV0 

e 


z(m); 


m=0 


Z  (m) 


(73) 
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and 


00 


A  00  .TT  1 

*  A  —  1  pyL.  ■ 

v~  ~L  (-1V)--  *<”>,  z(“°  -  (- l)me  4(f)  2  *_<*). 


m  — 


m=0 


(74) 


Here  (-iv)  ^  denotes  the  principal  value,  i.e.  (-iv)  ^  *  V  t 

We  now  multiply  the  incident  and  diffracted  waves  (73)  and  (74)  by 
an  arbitrary  function  a(v)  and  integrate  with  respect  to  V.  Then 
they  become 


and 


v  (t,  x)  ~  )  e  (0): 

o  Li  m 


(m) 


m=0 


*(t'  5>  ~  I  WW 


A  A  ( 


(m) 


m=0 


(75) 


(76) 


where 


er(t)  =  J  (-iv)-ra(v)e_i 


-ivt 


dV 


(77) 


We  note  that  (77)  holds  for  both  integral  and  fractional  values  of  r. 
Thus  (see  the  Appendix)  for  all  r 


e  (t)  =  I  e  (t) 
r  r  o 


(78) 


where  I  is  the  fractional  integral  operator  of  order  r. 
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By  comparing  (42)  and  (44)  with  (75)  and  (76)  we  conclude 


that  in  general 


««<*>  " 


(79) 


We  also  insert  (65)  into  (62)  and  compare  the  result  with  the 
A  (o) 

expression  for  z  obtained  by  inserting  (70)  into  (74).  Since 

z^  =  z  we  conclude  that 
o 

Do  -  •  (80) 

Summarizing  our  results  we  conclude  that  the  leading  term  of 
the  diffracted  progressing  wave  is  given  by 

V~e^(0)  (|-)W0)[t,  XjOi)]  | a ( 1  -t-  r %  ,  (81) 

where  0  =  c(t  -  T)  ,  p  is  given  by  (60)  and  D  =  D(T,  t),  a)  is  given 

A 

by  (71)  and  (67).  Thus  v(t,  x)  Is  given  parametrically  with  parameters 
(T,  T),  a)  by  (81)  and  (51). 
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SECTION  IV 


DIFFRACTION  BY  MOVING  EDGES 

The  progressing  wave  formalism  can  be  applied  to  problems  of 
scattering  by  moving  as  well  as  fixed  obstacles.  The  interesting 
problem  of  reflection  of  a  progressing  wave  by  a  moving  surface  is 
discussed  in  [2]#  There  it  was  found  that  the  usual  law  of  reflection 
is  modified,  and  a  Doppler  shift  of  the  instantaneous  frequency  is 
introduced.  It  is  also  interesting  to  note  that  if  a  time-independent 
progressing  wave  such  as  (5)  or  (43)  is  scattered  by  a  moving 
surface  the  reflected  wave  is,  in  general,  not  time-independent.  In 
problems  involving  moving  sources  [2,  3,  4]  general  (i.e.  not  time- 
independent.)  progressing  waves  also  arise  naturally. 

In  this  section  we  give  a  partial  treatment  of  diffraction  of 
a  progressing  wave  by  a  moving  edge  such  as  the  edge  of  a  finite 
cylinder  which  is  undergoing  translation  and  rotation.  As  in  the 
geometrical  theory  of  diffraction,  a  full  treatment  would  require 
the  solution  of  a  "canonical  problem",  involving  a  moving  infinite 
wedge.  However,  considerable  information  can  be  obtained  without 
solving  that  problem. 

Let  T|  be  an  arclength  parameter  on  the  edge  E  and  let 

X  -  £Cn,  t)  (82) 
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be  the  parametric  equation  of  E  at  time  t.  Then  is  the 

instantaneous  unit  tangent  vector  to  E.  We  introduce  the  vector 


s-it  ■  • 


(83) 


For  any  function  r|(t)  the  point 

x  =  x(t)  =  £[r|(t)  ,  t]  (84) 

is  constrained  to  remain  on  the  moving  edge.  Its  velocity  is 
•  • 

X  T|t^  +  5  and  the  normal  component  of  its  velocity  is 

£  -  +  it  -  Tltx  -  it-lili  ■  w*  (85) 

Therefore  we  refer  to  w  as  the  normal  velocity  vector  of  E  and  to 
its  magnitude  w  »  |w|  as  the  normal  speed. 

A 

As  in  Section  III  we  assume  that  the  phase  function  0  of  the 
diffracted  wave  (44)  is  equal  to  the  phase  function  0  of  the 
incident  wave  (42)  on  E,  Thus 

0Ct,  5_cn,  t>]  =  0[  t,  £  On,  t>]  ,  (86) 
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and  differentiation  with  respect  to  t  and  T)  yields 


A  A 

-  U)  +  k.5t  ■  -  0)  +  k.§t  on  E,  (87) 


A 

k.t^  *  k.t^  on  E  .  (88) 


If  we  multiply  (88)  by  -  and  add  the  result  to  (87)  we  obtain 


A  A 

-  UJ  +  k«w  *  -  U)  +  k.w  on  E 


(89) 


In  terms  of  the  unit  vectors  u 
become 


-  k  and  u  *  £  k  ,  (88)  and  (89) 
®  ID 


AA 

UJu-t^  *  U)u.  t^  on  E  ,  (90) 


—  u.w]  =  Ou[l 
c - 


—  u.w] 
c - 


on  E , 


(91) 


and  it  follows  that 


— 1 -  on  E 

-  1 
1  -  —u.w 
C - 


(92) 
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We  now  introduce  the  unit  vectors  N  «  —  w  ,  and  B  *  t , X  N  , 
Then  t^>  N,  B  form  an  orthonormal  system  but  N  and  B  are  not  the 

normal  and  binormal  vectors  in  general.  We  also  introduce  the  angles 

a  A 

3  ,  P  ,  €  ,  e  defined  by 


u  =  cos  0 

*1 

+ 

sin  0 

cos  € 

N 

+ 

sin 

e 

sin  e 

B  , 

(93) 

A  A 

A 

A 

A 

A 

u  =  cos  0 

*1 

+ 

sin  0 

cos  e 

N 

4- 

sin 

P 

sin  G 

B  . 

(94) 

Then  (92)  becomes 


A 

COS  0 


W .  a 

-  —  sm  p 


cos  0 


* 

COS  € 


1  -  —  sin  0  cos  € 
c 


(95) 


This  equation  may  be  called  the  law  of  diffraction  for  moving  edges. 

It  reduces  to  the  more  familiar  form  (49)  when  w  =  0,  In  general, 

A  A 

for  given  0,  e,  and  w  (95)  determines  0  as  a  function  of  C  (or  vice- 

A 

versa) ,  hence  a  one-parameter  family  of  directions  u  for  the  diffracted 
rays  emanating  from  E,  Furthermore (90)  and  (91)  yield 


where 


(96) 


cos  0 

A 

cos  0 


1  -  —  sin  0  cos  e 
_ c _ 

w  5  a* 

1  -  —  sin  0  cos  G 
c 


(97) 
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Since  U)  and  U)  are  proportional  to  the  "instantaneous  frequencies"  of 
the  incident  and  diffracted  waves  (see  [2,  3,  4]),  the  factor  ^ may 
be  viewed  as  determining  a  "Doppler  shift".  When  w  *  0  we  see  that 
1‘ 

It  is  interesting  to  note  that  back-scattering  from  the  edge 

A 

occurs  when  u  =  -  u.  In  this  case  we  see  from  (90)  and  (91)  that 


(u)  +  uu)  u._t  1  =  0 


on  E  , 


(98) 


and 


A  1  -  —  u.w 

n  c - 

U)  =  U)  - j - 

1  4*  —  U.w 
c - 


on  E  , 


(99) 


Since  w  =  |w|  <  c  we  see  from  (99)  that  uu  and  (JU  have  the  same  sign, 
hence  from  (98)  that 


_u .  =  0 


(100) 


Thus,  even  for  a  moving  edge,  back-scattering  occurs  if  and  only  if 
the  edge  is  orthogonal  to  the  incident  ray. 
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APPENDIX 


FRACTIONAL  INTEGRALS  AND  FOURIER  TRANSFORMS 

For  a  >  0  we  consider  a  class  C  of  functions  f(t)  which  vanish 
—  a 

for  t  less  than  some  real  number  t  and  for  which  |f(t)|^c  e^t  for 

o  1  o 

sufficiently  large  positive  t.  Then  the  Fourier  transform 

a(u>)  =  ~  J  eia,tf(t)dt  (101) 

is  analytic  in  I  uu  >  a  .  and 
y  m 

f(t)  =  J  a(UU)e"ia,tdO)  .  (102) 

Here  and  hereafter  the  path  of  integration  in  the  uu-plane  is  parallel 

to  the  real  axis  but  lies  in  the  region  I  uu  >  a  ^  0.  We  introduce 

m 

the  linear  operator 

fi  # 

1  f(t)  =  (-W))"m  a(0))e"ia,tdU)  ,  m  s  0  .  (103) 

m  v 

I  will  be  called  the  fractional  integral  of  order  m»  If  m  is  not 
an  integer  ( - iuu)  m  denotes  the  principal  value. 

Lemma  1 

t 

For  m  >  0,  I/Ct)  =  p^y  J  (t  -  T)m‘1f(T)dT  (104) 

-00 
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where 


00 


r(m)  *  (m-1)  1  * 

0 

Proof:  Let 

t  f  /  -iW'r 

I=J  (t  -  T)  e  dT 

-00 


/ 


-z  m-1 , 
e  z  dz 


(105) 


(106) 


and  let  UU  =  re^  where  r  =  |u)|  and  9  =  arg  U),  We  introduce  the  change 
of  variable  T  =  t  +  z/iUJ  in  (106).  Then 

„ei(e  -  TT/2) 

e_Zzm'1dz  (107) 


For  I  UU>0,0<9<tt  and  -  TT/2  <  9  -  TT/2  <  rr/2  .  For  values  of 
m 

0  -  TT/2  in  this  interval  it  is  easy  to  show  that  the  path  of 
integration  in  (107)  can  be  rotated  to  the  positive  real  axis  (i.e.  the 
integral  over  the  arc  at  infinity  vanishes).  Hence 


I  =  (-iUU)  me  ^*T(ra)  . 


(108) 


Thus 


r 

r(m)  J 


(t  -  T)m_1f(T)dT  1 


r(m) 


/ 


a(0U)  IdU) 


/ 


( - iuu) " ma (ou) e ” iaJtdou  =  I  f(t) 

m 


(109) 
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Lemma  2 


(110) 


n  m  n+m 

proofj  From  (103),  the  Fourier  transform  of  the  function  Imf(t)  is 
C-iUJ)  ma(uu).  Hence  from  (103) 

1  [i  f(t)]  =  I  ( - iuo)  n ( - itu)  ma(u))e  duu  =  I  f(t)  . 
n  m  J  rH-m 

(in) 
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